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The counterterm combination that describes the decay of pseudoscalar mesons into charged 
lepton pairs at lowest order in chiral perturbation theory is considered within the framework of 
QCD in the limit of a large number of colours Nc- When further restricted to the lowest meson 
dominance approximation to large- A^c QCD, our results agree well with the available experimental 
data. 
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1. The theoretical study of the tt" and rj decaying 
into lepton pairs and the comparison with the experi- 
mental rates offers an interesting possibility to test 
our understanding of the long-distance dynamics of the 
Standard Model. These processes are dominated by 
the exchange of two virtual photons and it is there- 
fore phenomenologically useful to consider the branch- 
ing ratios normalized to the two-photon branching ratio 
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with I3i{s) = -^/l — Am'j/s. The unknown dynamics is 
then contained in the amplitude A{Mp). To lowest order 
in the chiral expansion the contribution to this amplitude 
arises from the two graphs of Fig. 1 with the result 



A{s) = xp(m) + ^ 
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where = X») = -(xi + X2)/4 = X, with xi and X2 
the couplings of the two counterterms which describe the 
direct interactions of pseudoscalar mesons with lepton 
pairs to lowest order in the chiral expansion Q 

X [xMQbQrD^UU^ - QlQlD,,U^U) 

+ X2tr{U^QRD^UQL - UQlD^U^Qi,)] . (3) 

Here the unitary matrix U describes the meson fields and 
Ql = Qr = diagonal(2/3, — 1/3, — 1/3). The function 
C{s) in Eq. (|^) corresponds to a finite three-point loop 
integral which can be expressed in terms of the diloga- 
rithm function Li2(a:) — — jQ{dt/t) ln(l — t). For s < 0, 
its expression reads 
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FIG. 1. The lowest order contributions to the P £^£~ 
decay amplitude. The second graph denotes the contribution 
from the counterterm lagrangian of Eq. (^). 

The corresponding expression for s > Amf is obtained 
by analytic continuation, using the usual ie prescription. 
The loop diagram of Fig. 1 originates from the usual cou- 
pling of the light pseudoscalar mesons to a photon pair 
given by the well-known Wess-Zumino anomaly Q . The 
divergence associated with this diagram has been renor- 
malized within the MS minimal subtraction scheme of 
dimensional regularization. The logarithmic dependence 
on the renormalization scale fi displayed in the above ex- 
pression is compensated by the scale dependence of the 
combination xip) of renormalized low-energy constants 
defined above. Let us stress here that, as shown explic- 
itly in Eq. (|^) and in contrast with the usual situation in 
the purely mesonic sector, this scale dependence is not 
suppressed in the large-iVc limit, since it does not arise 
from meson loops. The evaluation of xit^) will be the 
main subject of this paper. 

It has recently been shown that, when evaluated 
within the chiral U{3)(S>U(3) framework and in the I/Nq 
expansion, the \AS\ = 1 transitions can also 

be described by the expressions (|l|) and (||), with an ef- 
fective constant Xk° containing an additional piece from 
the short-distance contributions Of course, a cast- 
iron understanding of these transitions is very important 
as the evaluation of x(m) could then have a poten- 
tial impact on possible constraints on physics beyond the 
Standard Model. We comment on this issue at the end 
of the paper. 

2. As a first step towards its subsequent evaluation 
we shall identify the coupling constant x in terms of a 
QCD correlation function. For that purpose, consider the 
matrix element of the light quark isovector pseudoscalar 
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density P^{x) = ^{uij^u — di"f^d){x) between leptonic 
states in the chiral limit. In the absence of weak interac- 
tions, and to lowest non-trivial order in the fine structure 
constant, this matrix element is given by the integral 

<i-{p')\p^{o)\r{p) > 
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with — ■|(2u7^w — djf^d ~ S7ps). In the chiral limit, 
the QCD three-point correlator appearing in this expres- 
sion is an order parameter of spontaneous chiral sym- 
metry breaking. This ensures that it has a smooth be- 
haviour at short distances. In particular, Bose symmetry 
and parity conservation of the strong interactions yield 
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withn{qlqliqi+q2f)=Hiqlqf,iqi+q2f). For very 
large (euclidian) momenta, the leading short-distance be- 
haviour of this correlation function is given by 
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Actually, what matters for the convergence of the integral 
in Eq. (^) is the leading short-distance singularity of the 
T— product of the two electromagnetic currents, which 
corresponds to 

hm n{{Xqf,{p'^p~Xq)\{p'-p)') 
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and which implies that the loop integral in Eq. (g) is 
indeed convergent. The QCD corrections of order 0{as) 
in Eqs. and ^ will not be considered here. Let us 
however notice that since the pseudoscalar density P^{x) 
and the single-flavour < il^ij} > condensate share the same 
anomalous dimension, the power-like fall-off displayed 
by Eqs. and is canonical, i.e. it is not modified by 
powers of logarithms of the momenta. 

On the other hand, at very low momentum transfers, 
the same correlator can be computed within Chiral Per- 
turbation Theory (ChPT). At lowest order, it is saturated 



by the pion-pole contribution, given by the anomalous 
coupling of a neutral pion, emitted by the pseudoscalar 
source P^(0), to the two electromagnetic currents, i.e. 
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where the ellipsis stands for higher orders in the low- 
momentum expansion and where Fq denotes the pion 
decay constant in the chiral limit. The matrix element 
< l~{p')\P^{Q)\l~{p) > itself may also be evaluated 
in ChPT. At lowest order, it is given by the diagrams 
of Fig. 1, where the (off-shell) pion is now emitted by 
the pseudoscalar source P^{0). The result reads, with 

t=(p'-p)2, 
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with the function A{t) defined in Eqs. (|) and (|). The 
contribution of the loop diagram of Fig. 1 is obtained 
upon replacing, in Eq. (||), the three-point QCD cor- 
relator by its lowest order chiral expression given in 
Eq. (H). The coupling constant x(/x) is thus given by 
the residue of the pole at i = of the matrix element 
< £~{p') I P^(0) I i~{p) >, after subtraction of the contri- 
bution of the two-photon loop, i.e. 
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Since the integral occurring in the above expression di- 
verges, we have regularized it by analytical continuation 
in the space-time dimension d. The coupling x(a*) on 
the left-hand side is then defined by the MS minimal 
subtraction prescription, as in Eq. (|^). Keeping only 
the contributions that do not vanish as {p' — p)^ goes 
to zero, and neglecting terms of order ©(m^/A^^), where 
Ah ~ Mp is a typical hadronic scale for non-Goldstone 
mesonic states, we obtain a somewhat simpler expression, 
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In order to perform this integral, one needs to extend 
the knowledge of the three-point correlation function 
"^{QitQ^Aqi + 92)^) in the chiral limit beyond its be- 
haviour at energies very high, Eq. (^, or at energies very 
low, Eq. (||). Stated like that, in full generality, this rep- 
resents a rather formidable task. As we shall next show, 
it is possible, however, following the examples discussed 
recently in Refs. to proceed further within the 

framework of the 1/iVc-expansion in QCD p^ . 

3. In the limit where the number of colours Nc be- 
comes infinite, with agXNc staying finite, the QCD spec- 
trum reduces to an infinite tower of zero-width mesonic 
resonances ||l4|, and the leading large-iVc contributions 
to the three-point correlator are given by the tree- 
level exchanges of these resonances in the various chan- 
nels, as shown in Fig. 2. This involves couplings of the 
resonances among themselves and to the external sources 
which, just like the masses of the resonances themselves, 
cannot be fixed in the absence of an explicit solution of 
QCD in the large-iVc limit. In this limit, however, the 
analytical structure of the three-point function in Eq. (||) 
is very simple: the singularities in each channel consist of 
a succession of simple poles. Furthermore, the quantity 
appearing in Eq. (|l2|) has the general structure 

lim ip'-pfniq^q^ip'-pf) = -'.^^ 
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where a priori the sum extends over the infinite spec- 
trum of vector resonances of QCD in the large-iVc limit. 
Equation (|l^) follows from the fact that its left-hand 
side enjoys some additional properties: i) In the pseu- 
doscalar channel, only the pion pole survives, while mas- 
sive pseudoscalar resonances cannot contribute, ii) The 
momentum transfer in the two vector channels is the 
same, iii) Its high-energy behaviour is fixed by Eq. (^. 
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Even though the constants ay and by depend on the 
masses and couplings of the vector resonances in an un- 
known manner, they are however constrained by the two 
conditions 
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which follow from Eqs. (|9|) and (||), respectively. No- 
tice that there are no contributions from the perturba- 
tive QCD continuum to these sums. Taking the first of 
these conditions (which, coming from the anomaly, has 
no 0{as) corrections) into account, we obtain 
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This equation, together with the two conditions (|l4|), 
constitutes the central result of our paper. This is as 
far as the large-iVc limit allows us to go. Let us point 
out that the scale dependence of x(/i) is correctly repro- 
duced by the expression ( |l5|) , again as a consequence of 
the first condition in Eq. (|l4|) . However, in order to ob- 
tain a numerical estimate of xip-) additional assumptions 
are needed. 

4. In order to proceed further, we shall consider the 
Lowest Meson Dominance (LMD) approximation to the 
large- A^c spectrum of vector meson resonances discussed 
in Jisf . This approximation has been shown to reproduce 
very well the relevant low-energy constants of the 0{p'^) 
chiral Lagrangian |1^] and the electromagnetic 7r+ — vr" 
mass difference |ll]. In our case, it corresponds to the 
assumption that the sums occurring in Eqs. (|l^) and ( p^ ) 
are saturated by the lowest lying vector meson octet. In 
the LMD approximation to large-iVc QCD, the two con- 
ditions ( p^ ) completely pin down the two quantities ay 
and by in terms of Fq and of the mass My of this lowest 
lying vector meson octet. 
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In fact, within the LMD approximation of large- Ac 
QCD, it is easy to write down the expression for the 
correlation function 7i(g^, g|, {qi -\- 92)^) which correctly 
interpolates between the high energy behaviour in Eq. (|^) 
and the ChPT result in Eq. (|) @ 
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FIG. 2. The contributions to the vector-vector-pseudosca- 
lar three-point function in the large- A^c limit of QCD. The 
sum extends over the infinite number of zero-width vector (V^) 
and pseudoscalar (P) states. 
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Notice that this expression also correctly reproduces the 
behaviour in Eq. (||). With the results of Eq. (|6|), and 
for Nc='i, it follows from Eq. (|l|) that 

LMD( ) = 11 _ 3 / 4 2il? ^ (18) 
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Numerically, using the physical values Fq 
and My = M„ 



p - 770 MeV, we obtain 

= 2.2 ±0.9, 



92.4 MeV In view of these uncertainties we conclude that it does 
not seem to be possible, within our understanding of 
long-distance effects in the electroweak interactions, to 

(19) 

argue that M^M is, at present, a useful decay to 



where we have allowed for a systematic theoretical error 
of 40%, as a rule of thumb estimate of the uncertainties 
attached to the large- TVc and LMD approximations. The 
predicted ratios of branching ratios in Eq. (|^) which fol- 
low from this result are displayed in Table |[ We 
conclude that, within errors, the LMD-approximation 
to large-iVc QCD reproduces well the observed rates of 
pseudoscalar mesons decaying into lepton pairs. 

5. At present, the most accurate experimental deter- 
mination of the K1 l^^l-i-^ branching ratio [20] gives 
the result: Br{Kl ^1+^-) = (7.18 ± 0.17) x 10"^. 
In the framework of the 1/Nc expansion and using 
the experimental branching ratio [2] Br{K^ 77) = 
(5.92 ± 0.15) X 10~^, this leads to a unique solution for 
an effective Xk'I — 5.17 ± 1.13. Furthermore, following 
Rcf. [5], XkI = X - a/" Sxsd where TV = (S.e/gsCrod) 
normalizes the —> 77 amplitude. The coupling gg, 
governs the AI —1/2 rule, the constant Cred is defined in 
Ref. [5] and SxI'd""""' = (+1.8 ±0.6) is the short distance 
contribution in the Standard Model [6]. 

Therefore our understanding of the short distance con- 
tribution to this process completely hinges on our under- 
standing of the long distance constant J\f and therefore 
of the AI =1/2 rule in the 1/Nc expansion. Moreover, 
Cicd is regretfully very unstable in the chiral and large- 
Nc limits, a behaviour that surely points towards the 
need to have higher order corrections under control. For 
instance, for A/^r = 0, Mk ^ 0, Nc 00 one obtains 
Crcd — 0, while for Af^ = Mk — 0, Nc 00 (and the 
external off shell) one obtains Cj-cd — —4/3 instead. 
The analysis of Ref. [5] uses Crcd — +1 and gs — 3.6, 
where these numbers are obtained phenomenologically 
by requiring agreement with the two-photon decay of 
Ar°,7r°,77 and rj' as well as K ^ 27r,37r. Should we use 
these values of Cred and gg and Eq. (19) we would ob- 
tain XkI ~ ^'^ corresponding to a ratio R{K^ — + 
fi+^i-) = (2.24 ±0.41) X 10"^ which is 2.5cr above the ex- 
perimental value RiKl /i+/i") = (1.21 ±0.04) x 10"^ 

TABLE I. The values for the ratios R{P £+£-) obtained 
within the LMD approximation to large-Ac QCD and the 
comparison with available experimental results. 
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constrain physics beyond the Standard Model. 
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